wo: 



Remarks on generalized scalar-tensor models of 

dark energy 

M. Alimohammadi* and H. Behnamian 

Department of Physics, University of Tehran, 
" North Karegar Avenue, Tehran, Iran. 

o 
o 

(N 

_ ■ Abstract 
Qh. 

D , The generalized scalar-tensor models with Lagrangian F((f>, R)—U(cj))(\74>) 2 

^0 ' are considered. It is shown that the phantom-divide-line crossing and the 

| deceleration to acceleration transition generally occurr in these models. 

■ Two specific examples, the coupled quintessence model and the Brans- 

Dicke model are considered. For the first example, it is shown that for 
the models with £ > 3/16, the ui = —1 transition exists. This is verified 
numerically for some special cases. For the Brans-Dicke model, it is shown 
that the transition does not occur, a result which can be verified by using 
the exact solution of this model. Finally the contribution of quantum ef- 
fects on these phenomena is investigated. It is shown that for some special 
cases where the u) = — 1 transition is classically forbidden, the quantum 
£\J ■ effects can induce transition. The £ = 1/6 of coupled quintessence model 

^ ' is an example of this. The quantum effects are described via the account 

C , of conformal anomaly. 

; 1 Introduction 

Nowadays based on various observational data, it has been verified that our 
universe is now in an accelerating phase [1]. The source of this accelerated 
expansion, which is known as dark energy, is believed to compose nearly 70% 
of the present universe. The precise nature of dark energy is not yet known 
but all of the models describing the dark energy have a common property: they 
produce the repulsive force, or in other words, they have the negative pressure. 
Introducing the equation of state parameter by u> — p/p, where p is the pressure 
and p is the energy density, it must satisfy u> < —1/3, in Priedmann- Robertson- 
Walker (FEW) background metric, to ensure the accelerated expansion. 

One of the interesting features of dark energy is the dynamical behavior of 
equation of state parameter lo. This is because some astrophysical data slightly 
favor an evolving dark energy and show a recent to = — 1, the so-called phantom- 
divide-line, crossing [2], To describe the behavior of u>, many dynamical dark 
energy models have been introduced. The simplest of these models are the single 
component scalar models, including the quintessence model, which consists of 
a normal scalar field [3] , and the phantom model, which is a scalar field theory 
with unusual negative kinetic energy [4] . None of these models can describe the 
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u> = — 1 crossing. In quintessence model u> is always u> > — 1, and in phantom 
model cj always satisfies u> < — 1. A possible way to overcome this problem 
is to consider the quintom model [5], the model consists of one quintessence 
and one phantom fields. It can be shown that the transition from u> > — 1 to 
u> < — 1 always occurs in the quintom model with slowly- varying potential [6] . 
The interaction of single-component scalar field with background (dark) matter 
can also induce this transition [7] . 

An alternative candidate of dark energy are the models known as modified 
gravities, in which it is postulated that the gravity is being nowadays modi- 
fied by some extra terms which grow when the curvature decreases. The f(R) 
gravity, whose action is a general function f(R) of the Ricci scalar R, is the 
most famous modified gravity theory [8]. The /(G) gravity and its general- 
ization F(R, G) gravity, where G being the Gauss-Bonnet invariant, are also 
introduced [9]. Several aspects of these models, including their behavior un- 
der lo = — 1 transition, have been studied recently. It has been shown that 
the phantom-divide-line crossing and the deceleration to acceleration transition 
generally occur in F(R, G) models. The contribution of quantum effects on 
these phenomenon has been also obtained [10]. 

An important class of modified gravity theories is the scalar-tensor theories. 
In these models , there exists a scalar field 4>, besides the usual space-time metric 
<7 M „, to describe the gravitational interaction. The earliest type of these models 
is the Jordan-Brans-Dicke theory with the following action, in the so-called 
Jordan frame, [11] 



Here we use the units in which h = c = 8irG — 1. 'w' is a dimensionless 
parameter and S m is the action of dustlike matter. The scalar field <p plays the 
role of the inverse of the gravitational coupling, which is constant in general 
relativity but in Jordan-Brans-Dicke theory is variable. For the case where 
V{4>) — 0, the model is called the Brans-Dicke theory. 

The generalization of the Jordan-Brans-Dicke theory is the scalar-tensor 
theory which is described, in the Jordan frame, by the action [12] 



Several aspects of scalar-tensor models have been studied, including their be- 
havior as a dynamical system [13], the gravitational wave and inflation in this 
framework [14], and their perturbative aspects [15, 16]. Two review articles on 
this subject are [17]. 

The most general action in the Jordan frame, based on the scalar-tensor 
theory, is the following: 



which is called the generalized scalar-tensor (ST) theories [16,17]. The scalar 
theories, i.e. quintessence and phantom models, f{R) gravity models, Jordan- 
Brans-Dicke theory ([1]), and the scalar- tensor models © are special examples of 
action ([3]) . It is must be noted that the terminology is not fixed in scalar-tensor 
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theories. The term "generalized scalar-tensor models" has been used for less 
general actions than eq.©, even for action with F(<f>) = <f>. 

The present paper is devoted to the study of generalized ST theories. We 
want to study the conditions under which the lo = — 1 crossing occurs in gener- 
alized scalar-tensor theories. The condition of transition between deceleration 
and acceleration phases is also studied. We show that, through some specific 
examples, the phantom-divide-line crossing can not occur for some models. It 
is shown that by considering the quantum contributions, the situation becomes 
different. In fact, some of the forbidden transitions become allowed, i.e. the 
lo = — 1 transitions, for these cases, are quantum induced. The same phe- 
nomenon has been seen for quintessence and phantom theories [18]. The quan- 
tum effects arc described via the account of conformal anomaly, reminding about 
anomaly-driven inflation [19]. The contribution of the conformal anomaly in en- 
ergy conditions and Big Rip of phantom models has been discussed in [20] . 

The scheme of the paper is as follows. In section 2, the Friedmann equations 
of generalized ST theories for FRW space-time are obtained. By perturbative 
solving of these equations around the to = —\ and a = transition points in 
sections 3 and 4, respectively, the conditions of occurrence of these crossings are 
obtained. In section 5, we consider the explicit example of coupled quintessence 
model in which there exists a free parameter £. It is shown that the lo = — 1 
transition occurs if £ > 3/16. It is seen that the numerical solving of Friedmann 
equations verifies our perturbative results. In section 6, the Brans-Dicke model 
is considered as a second example. Our method shows no transition for this 
model, a result which can be verified by exact solution of this model. Finally 
in section 7, we introduce the quantum contributions (comes from conformal 
anomaly) to this problem, and show that it can induce the to = — 1 transition 
for special cases which are classically forbidden. For example for £ = 1/6 case of 
the coupled quintessence model, in which the transition is classically forbidden, 
the quantum phenomena induce the lo = — 1 transition. 



2 The generalized scalar-tensor models 

Consider the generalized scalar-tensor gravity with action ©. Varying this 
action with respect to the metric g^ results in [17] 

(iV - V„V„ + g^U)F R (4>, R) - hg^ [F(4>, R) 

-u{4>){V4>) 2 } - u{4>)v^4>v v 4> = T^. 

In this equation, T M „ is the energy-momentum tensor of matter field 

2 5S, 

and 

F R {<f>,R) = — ^ — , F <t ,(<j>,R) = — — — . (6) 
Varying the action with respect to the scalar field tfi, results in 

F+ty, R) + 2tf(0p0 + 2VVV M tf(0) - ^W) 2 - 0. (7) 

a<p 



1^ = ^=^, (5) 



3 



For the case where U{<f>) = rj =constant, eq.([7]) reduces to the corresponding 
equation of motion in [17]. 

A spatially flat FRW space-time in comoving coordinates (t, x, y, z) is defined 
through 

ds 2 = -dt 2 + a 2 (t){dx 2 + dy 2 + dz 2 ), (8) 

where a(t) is the scale factor. For this metric, the (t, t) component of the 
evolution equation ((¥]) for homogenous scalar field <p(t) becomes: 

'SH 2 F R {4>, R) = Pm + \u{^)4> 2 - 3HF RR (<t>, R)R 

(9) 

- 3HF R(i> (^ R)j> - ±F(<f>, R) + ^RF R {ct>, R). 

H(t) = a(t)/a(t) is the Hubble parameter and p m is the matter energy density 
with the evolution equation 

p m +3H(p m +p m ) = 0. (10) 

R is expressed in terms of H as following 

R = 6{H + 2H 2 ). (11) 

The (x,x) component of eq.(j3]), when eq.® is subtracted from it, becomes 

-2HF R {<I>, R) + HF RR ((f>, R)R + HF R<P {^ R)j> - F RRR {<j), R)R 2 - 2F RR4 {<I>, R)R<j> 

- F RR {cf>, R)R - Fr^, i?)0 2 - F m (<f>, R)4> - £/(0)0 2 = lmPm . 

(12) 

In above equation, j m is defined by j m = 1 + u) m , where uj m = p m / p m , and the 
subscripts denote differentiations, e.g. 

F RRt j,(4>,R) = d2R94> . (13) 

Finally the scalar field equation of motion (J7|) in FRW metric becomes: 

Ff(<p, R) - 2U(<f>)(4> + 3H<j>) - f/(0)0 = 0. (14) 

In above equation U = (dU/d(j))4>. Equations (0, fTU]). (fHf and ([11]) are the 
Friedmann equations of generalized scalar-tensor models. Note that these four 
equations are not independent, in fact, taking the time derivative of eq.([5]), 
adding it cq. lfl^]) . multiplied by 3H, and using eq. tfTI"]) . result in cq. lflU]) . So 
there are three independent Friedmann equations. In the case of ordinary 
scalar-tensor models where F(<fi,R) = RF(<fi) — V(<p), the resulting Friedmann 
equations lead to the corresponding ones in scalar-tensor theories. 



3 The phantom-divide-line crossing 

For the ordinary dark energy models and in the context of Einstein gravity, the 
equation of state parameter uj = p/p is obtained via u> = —1 — 2H/3H 2 . For 
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other theories, including the usual and generalized ST theories, the effective 
equation of state parameter Lu c g is also defined through [9, 16] 



2 H 

W e g = -l-g^2. (15) 

So u e g crosses the phantom-divide-line, i.e. ui e s = —1, if H(t) has a relative 
extremum at some time t = to- Restricting ourselves to t — to -C h^ 1 , where 
ho = H{to) and h^ 1 is of order of the age of our universe, the Hubble parameter 
can be expanded as 

Hit) =h + h 1 (t- to) a + h 2 (t - t Q ) a+1 + 0((t - t ) a+2 ), (16) 

in which a > 2 is the order of the first nonvanishing derivative of H(t) at t = to, 
and hx = (l/a!)iJ( Q )(t ). ffW^) is nth derivative of H(t) at t = t . The 
transition occurs from lu > — 1 to < — 1 regions when a is an even positive 
integer and hi > 0. For hi < 0, and even integer a, the system goes from 
ui < — 1 to > — 1 at t = to- We now examine whether there is any solution 
for Friedmann eqs.(9), (10), and (14), if H(t) is given by (16). If for an even 
positive integer a, we can find a nonzero value for hi, then the possibility of 
u> = — 1 transition of generalized ST theories is proved. 
We first consider the Friedmann equation ((9]): 

H 2 F R (^R)=p(t), (17) 

in which 

lr 1. 



= §[ Pm + ^U{4>)tf - 3HF RR (<t>, R)R 
SHFtytt, R)4> - \F{cp, R) + ^RF R (<p, R) }. 



(18) 



Expanding both sides of eq. ([T7j) around to = 0, one finds, up to order t 2 ,: 

h 2 F R (0) = 0(0), (19) 
h 2 (F RR R + F R4 ^) t= o = 0(0), (20) 



and 



-jho(F RRR R 2 + 2F RR< j,R(j) + F RR R + Fn^cf) 2 4- F R ^4>)t = o 
+ 2/io/iiF H (0)* a ,2 - -0(0). 



(21) 



The second Friedmann equation ()T0|) . at zero order, leads to 

Pm(0) + 3/lo7mPm(0) =0, (22) 
and the last equation of motion (|14p . up to order t, results in: 

F (O) - 2C/(O)(0 + 3M)t=o - C/ (O)0 2 (O) = 0, (23) 

and 

^(0)0(0) - C7*(0)(4# + 6/i o 2 ) t= o - 2t/(O)(0 + 3M)t=o 
- C/ 00 (O)0 3 (O) + 12hiF m (0)8 a , 2 = 0. 



(24) 
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Using egs.flTTj), (EU) and , the eqs.QU) and (gDJ) lead to: 



3h 2 F R (Q) + Pm (0) + it/(O)0 2 (O) - 1f(0) 
- 3/10^(0)^(0) - 36^1^(0)^,2 = 0, 

and 



6[3/i (/i2 + h hi)F RR + \2h h\F RRR + 2h h 1 F RR ^] t=0 8 aa + 18h hiF RR (0)S a . 3 

(26) 



~[ 7mPm - h Q F R(t> (f) + U{cj))4> 2 + F R4>4 ,<\) 2 + F R( p(f> ]t= ho = 0, 



respectively. Note that at i = 0, -R = R(0) = 12/i§. 

Let us first consider the case a = 2. The eqs. (|2"2"j) . (|25p and (J2SJ) then result 
in /io, ^i and /12, respectively, ho and /ii are 

*o = -J*^ t (27) 



37mPm(0)' 



and 



fel = ifih I 7rt\ Sh ° FR + pm + ^W^ 2 - ^ F - 3^0^.00 ]t=o. (28) 

So the parameter /ii is, in general, different from zero, which depends on initial 
values p m (0), p m (Q), 0(0), and 0(0), can be positive or negative. This proves 
that, in general, the generalized ST theories can explain the phantom-divide- 
line crossing. The remaining two equations (|2"3")l and (f2"4"|) determine the values 
of 0(0) and 0(0) in terms of the four mentioned initial values. The higher order 
terms of the expansions specify the other coefficients /i3,/i4,... For example 
eq. (|2~T]) gives the parameter h 3 . 

It is clear that when the Lagrangian F((j>, R) is linear in R, as is the case 
for scalar-tensor theories with action then eq. (|2"8|) is not the solution (note 
that in this case, the last term of eg. (j2"S"|) vanishes). For these cases, it can be 
shown that eg. (|2"Tj) specifics hi as follows: 

hl = a J~ \n\ fiPmlfnho + 6h U{(j))(j) 2 - F + F R4> (h 4> - <j>) 

4i< Jj(0) (29) 

+ F R(j>( j,{h Q ^ 2 - 300) - -FR0000 3 ] t= o- 



For the next choice a — 3, it is clear from eg. ([25]) that it does not lead to 
an expression for hi, unlike the case a = 2 in which eg. ([25]) leads to (|28|) . but 
instead, this eguation leads to an extra relation between the initial values p m (0), 
p m (0), 0(0) and 0(0). For a = 3, eq. ([2"6"|) determines hi. So the a = 3 solution 
is also possible whenever there exists a special relation between the four initial 
values. The situation is worse for a > 4, since in these cases eg. (|2l)f also leads 
to another relation between the initial values. So except for these fine-tuned 
initial values, the only acceptable solution is a = 2, which always permits the 
u) = — 1 transition. 
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4 The deceleration to acceleration transition 



To study the transition from a < to a > 0, we first note that 



- = H + H\ 

a 



(30) 



So at time t' = where a = 0, one has £T = — H 2 . Expanding H(t) around 
t' = 

H{t) = H Q + Hxt + H 2 t 2 + ■■■ (31) 



then H = —H 2 



results in 



t=o 



Hi 



-Hi 



(32) 



Similar to the previous section, one can seek the solutions for eqs.©, (fit)]) and 
(fT4|) . when H(t) is given by (f3Tjl . with condition (|32| . In this case, the resulting 
relations for Hq and i?2 are 



Hn 



Pm(0) 
37mPm(0) 



(33) 



and 



36^0^(0) 



i[/(0)0 2 + 72i7 4 ^ - 3H F R ^ - 



(34) 

respectively. The other coefficients of expansion (|3"Tj) can be also found consis- 
tently, which proves the existence of a — crossing in generalized ST theories. 
For the cases where F(4>, R) is linear in R, H 2 in eq. ([3"4")) is not more the solution 
and another expression, like eq. (|29jl for uj = — 1 crossing, can be found. 



5 The coupled quintessence model 

As an explicit example, we consider the coupled quintessence model introduced 
by following action: [17] 

S = \ J d*x^d [(1 ~ t4?)R - (V0) 2 ] + S m , (35) 

in which we do not consider the potential V(<p). At ^ = 0, cq. (l33|) becomes 
the ordinary quintessence model, and for £ 0, the scalar field <fi coupled non- 
minimally to the gravity. The reason for considering the parameter £, first 
introduced in [21], has several answers which have been discussed in [17]. Per- 
haps the main reason to include a £ ^ term is that it is introduced by the 
first loop corrections and is required by normalizability of the theory [22]. If 
one considers a classical theory with £ = 0, renormalization shifts it to one with 
£ 0. We will return to the quantum corrections in section 7. 

In this section, we study the conditions under which the lo = — 1 transition 
occurs in V((j)) = coupled quintessence model. For a — 2, our zero-order 
relations are eqs. ((2"2"|) , (j2"3"|) and (|25[) . Since here we will not consider the matter 
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6H 2 (1 - £,cj) 2 ) - 2 - UH$<M> - 2p = 0, 
- 2H(1 - t4> 2 ) - 2HtM> + 2^4> + (2£ - 1)4> 2 =p + p, 



field, for simplicity, the eq. ([2"2")) . which is the Friedmann equation (TlT)]) at zero 
order, is not suitable. Instead we consider eq. (fT2"|) and our three independent 
Friedmann equations become eqs.Q, (JT2J) and ([U]), where for action (j3"5|) . result 
in: 

(36) 
(37) 

(38) 

(39) 

(40) 
(41) 



and 



4> + 3i?0 + £<p(6H + 12H 2 ) = 0, 
respectively. At zero order, these equations lead to: 

<j> 2 + 12^o0o0o - 6^(1 - C0o) = 0, 

- 2£/i o 0o0o + ^<j>o4>o + (2£ - 1)4>1 = 0, 
0o + 3/i o 0o + 12^o^o = 0, 



respectively. Here we assume that S m = 0. The first equation specifies 4>q 
in terms of (fio, and by eliminating 0o between eqs. (|40[) and (|41|) . one more 
expression for 0o, in terms of (f>o, is obtained. These relations are 



-h 



and 



6^o±v 36 ^o+6(i-e0i) 
/'o 



»o = 



l-2£ 



4^o ±^0^8(6^-1) 



(42) 



(43) 



respectively. Note that eq. (|4"3"|) . which comes from eqs. (j3"T|) and (|38l) . holds only 
at the uj = — 1 transition point. This is because by setting H = 0, which is 
the signature of to = —1 transition, eqs. ([37| and (|38ll lead to (|40| and (j4"Tjl and 
then (|4"3"|) is found. This is not the case for eq. (|4"2")) which holds at any instant 
of time. 

The reality constraint of the scalar field </>, restricts £ in cq. (|4"3"|) by 



<4 

Setting equal the eqs. (j4*2"|) and (j4*3"j) . results in 



(44) 



-1 [ v^-li^J =±(1-2^^36^ + 6(1-^) • (45) 

Note that for £ > 1/2, the left-hand-side of eq. (l4"5|) is positive, so in the right- 
hand-side, only the minus sign must be used for £ > 1/2. Solving (f>o from 
cq. (|45|) . one finds 



, (±) 2 _ 3(1 -2Q 2 



2£ 2 



(6£ - 1) (6£ + 1 ± 2^12^-2) - 9(1 - 2£) 2 



3(1 - 2Q 2 
2£ 

(46) 

(-)' z 

It can be easily seen that <fi is always negative for all £ > 1/6, so there 

(+) 2 (-M 2 
9 , it can be seen that O < for 



n -1 



exists no real solution for 



For 
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1/6 < £ < 3/16 and 0g > for £ > 3/16. So the real value of the scalar field 
(j) at transition time t = 0, is from eq.(|46p. and it exists if 

« > ^ ^ 

In other words, the coupled quintessence model can cross the u) = — 1 line only 
when £ satisfies eg. (|4"T| . This leads us to consider the eq. (f4"2"|) and the right- 
hand-side of eq.gSJ) with "+" and "-" sign, for 3/16 < £ < 1/2 and £ > 1/2, 
respectively!]]. 

To determine the crossing behavior of this models, we must calculate the 
parameter hi from eq. (|29p . or instead, from perturbative calculation of eq. (|12[) . 
or (pT7|) . To do so, it is better to rewrite the differential eqs. ([3"6"|) - ((3"5|) in terms 
of red-shift parameter z, defined through 

l + z = —, (48) 
a 



instead of the time variable t. In eq. (|48j) . a* is the scale factor of the universe 
in the present time. Using 



d 2 2 , x2 d 2 , , 2 dH d o, . d <4!,) 
— = H 2 (l + z) 2 — + H(l + z) 2 — — + i/ 2 (l + z)—, 
dr z dz z dz az az 



the first Friedmann equation (|36p results in 77(2) = d^>/dz in terms of 4>(z), and 
eqs.([5T|) and specify dH/dz and d^/dz. The resulting equations are: 

„(*) = — , (50) 



(1 + z)ri(z) = 6£0(z) ± y/36e^ 2 (z) + 6(l-^ 2 (z)) , (51) 



dff(z) _ 24£V(z) + (1 - 2£)(1 + z)VM - 8£(1 + z )rj{z)^z)_ ^ 



2(l + z)[l + (6£ 2 -O0 2 W 

x {24£ 2 3 (z) + (2f - 1)(1 + z)V(z) 



dz 2(l + z)[l + (6£ 2 -£)0 2 (z) 
+4(1 + z)[ 1 - (£ + 12£ 2 )0 2 (z) ]t7(2) + 2(1 + z) 2 (7£ - 6$V(*M*) - 24£0(z)} 

(53) 

As explained after eq. (|4"Tl) . for 3/16 < £; < 1/2, cq. (l5Tj) must be used with " + " 
sign and for £ > 1/2, the " — " sign is accepted. 

Expanding H(z) around the transition point zq, one finds 

H(z) = h Q + h' 1 (z~z ) 2 + h' 2 (z-z ) 3 + --- . (54) 



1 Since only the plus sign of eg. 1146 \ , or eq. ((43} and the LHS of eq. (J45J , leads to acceptable 
solution for <fio, we must consider the eg. 1143 II and the LHS of eq. {(45} with " + " sign. Under 
this condition, the LHS of eq. (1451 1 becomes positive for 3/16 < £ < 1/2, therefore the " + " 
sign must be chosen in the RHS of this equation and therefore in eq. l|42jl . For § > 1/2, the 
accepted sign in the RHS of eq. (I45P and in eo,l !42H is " — ", as pointed out after eq. )45t . 
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in which we choose a = 2. Note that oj = -1~{2/3)H/H 2 = -1 + ^^-dH/dz, 
so at transition point zq, one has AH / 'dz\ z=Zo =0. It is clear from eg. (|52() that 
dH/dz = A(z)H(z), so A(z ) = from which d 2 H/dz 2 \ z=Zo = H(z)dA/dz\ z=ZQ . 
In this way, eas.([5l3 ]) -(|53" j) finally result in h[ = (l/2)d 2 H/dz 2 \ z=Zo as follows 



ho 



x {96(i + zo)[ (e^ 4 +e+ e)4>i + ?e ~ e 



o V l + ^) 2 [l + (6e 2 -0^] 2 ^ L ^'^ lV " S ,S ' Sm) 
4(1 + 2 ) 2 [ (36£ 3 - 36£ 2 - 30^»o - 10£ + 3 } V 2 + 12(1 + z ) 3 (4£ 3 
2(1 + z ) 4 (4£ 2 - 4£ + 1) % 4 + 192£ 2 2 - 192£ 3 % 4 }. 



6£ 2 



- 3^0^ 
(55) 



In this equation, <po — 4>(zq) and rjo = 77(20)- 

At the first step, we solve numerically the differential equations (j5U|) , ([52")) and 
([53"]) for the specific example £ = 1. By choosing the initial values <p(z = 0) = 1/2 
and H(z = 0) = 72 and calculating 77(2 = 0) from eq. lfST]) with " — " sign, because 
£ > 1/2, one finds H{z), 4>{z), rj(z) andw(z) = -1 + ^p-dH/dz in Figs.(l)-(4). 
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Figure 1: The plot of H(z) of coupled quintessence model with £ = 1. 
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Figure 2: The plot of <f>(z) of coupled quintessence model with £ = 1. 
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Figure 3: The plot of rj(z) = d<fi(z)/dz of coupled quintessence model with 
£ = 1. 
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Figure 4: The plot of ui(z) of coupled quintessence model with £ = 1. It crosses 
uj = — 1 at z = 0.393. 



As it is clear from Fig. (4), uo(z) crosses the u> = — 1 line, from u> > — 1 region 
to uj < — 1 region, at z — 0.393. At this point, H(z ) = 112.83, (j>(zo) = 0.16, 
and t)(zo) — —1.18. 

Our perturbative results verify this crossing. For £ = 1, ^ +) from eq.® be- 
comes </>q + ' > = 0.16, which is in coincidence with numerical calculations. Choos- 
ing zq = 0.393 and ho = 112.83, eq. (p)Tj) results in r/(zo) = —1.18, and h[ from 
eq. (f55|) becomes h{ — —618.9. Therefore, up to lowest order, we have: 

2(l + z)dH Ah\ , 

3H{z) dz 3h (56) 

= -1 - 7.31(1 + z){z - 0.393) + • • • 

Fig. (5) shows the relation (|56"|) . which has the similar behavior, near the transi- 
tion point z — zq, as the Fig. (4). Note that in both approaches, the number of 
needed initial values are the same. In numerical calculation, we use two initial 
values 4>{z = 0) and H{z — 0), and in perturbative calculation, two parameters 
zq and ho = H(z — zq) are chosen. 
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w 




-3 - 



Figure 5: The perturbative variation of lu(z) of coupled quintessence model 
with £ = 1. 

6 The Brans-Dicke model 



As a second example, we consider the Jordan-Brans-Dicke model with V((j>) — 
and S m = 0, 



S= \ /<!', , 



(57) 



Comparing (|57| with (J3]), gives F{(j)) = Rcj) and U(4>) — w/</>. In this way, 
cqs.([23])-([M]), with a = 2, and (J2S]) become: 



+ wA 2 - 2w(3h A + B) = 0, 

C - A 3 + 2AB + 3h A 2 - 3h B + —hi = 0, 

w 

- wi 2 + 6h Q A + 6hl - 2— = 0, 

00 

B + wA 2 -ftoA+^±£°=0, 



(58) 
(59) 
(60) 

(61) 



and 



hi 



-C + h B + 6h wA 2 - Ylh\A 



3 Pol 2 h 



respectively. In above equations, the parameters A, B, and C are 



(62) 



A= Z 



B= ^ 



C= ^ 



(63) 



and R(0) — \2h\ has been used. Setting p = p — in eqs. 
S m — 0, eq. ([6*0|) then gives A as follows: 



to ensure 



3h 



w 



2w 



Reality of A demands 



w > 



2' 



(64) 



(65) 
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where for the reasons which will be discussed later, we do not consider the 
w = —3/2 case. 

The parameter B can be found from either eqs. (j55|) or ([BTj) . Using eci. ([l)U|) . 
they yield 

B = -(5h A + %h%) = —hi + 3(- - l)h A. (66) 
w w 

Putting dni| into eq. ((55|) . the parameter V is fixed as 

w = -|. (67) 

In other words, the Brans-Dicke theory has lo = —1 only when w = —4/3. 
Putting ([67]) back into eq. (l64|) . with plus sigrH, and eq. ((66|) . results in 

A = -3/i , B = 9/iq ■ (68) 

The parameter hi can be found by either eqs. ([I)9"|) and ([6"2"|l . where using |6T|) 
and (|68p. it gives the result 

fci = -|(C + 27Ag) = -i(C + 27Ag). (69) 
The above equality gives C: 

C = -27hl, (70) 

from which 

fti = 0. (71) 

So it seems that there is no tu — — 1 crossing phenomenon in Brans-Dicke model. 
In fact, calculating the other coefficients of expansion (fT6|) results in hi(i > 1) = 
0. 

To justify our result, let us exactly calculate the equation of state parameter 
u> of the Brans-Dicke model. We first write eqs.©, (fT2")l and (fT4^1 for this model, 
which results in 

3# 2 </> = - 3#0, (72) 

- 2H<t> + H<j>- <j>- -cj> 2 = 0, (73) 

and 

n W " * W ' o 

6( J ff + 2i/ 2 )-2-(0 + 3 J ff0) + — 2 = O, (74) 

<p <r 

respectively. Using the red-shift parameter z, the above equations can be used 
to obtain the following relations in z-space: 

6 ="< i+ ' )2 6sy +6 < i+i >if <75) 



dH{z) H(z) 



dz (l + z)(6 + 4w) 



w( 1 - - 2w H I + z) 2 ( ^tt] + 8w(l + z)^4 + 12 
<P(z)/ <P(z) 

(76) 



2 Note that for the minus sign of eg. 1164 \ , the solution of eg. (l66l > is w = —3/2, and for the 
plus sign it is w = —3/2 and w = —4/3. Because of the physical condition l|65| l. w = —4/3 is 
the only acceptable solution and therefore we only consider the plus sign of eg . J61H . 
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and 



drj(z) 4>( z ) 
dz ~ (l + z)(6 + 4w) 



-w(l + 2w)(l + z) 



r,*{z) 



(77) 

-12w(l + 2 )^ + 8(w-3)^ + -^ 
<p z {z) <p{z) 1 + z 

in which r/(z) — d<j>(z)/dz. Because of the denominator of cqs. ([7S)) and ([77| . 
we take w ^ —3/2 for these equations. We will consider the w = —3/2 case 
separately. Solving eq. (|75p for rj(z)/(j)(z), results in 

n{z) i 



(z) w(l + z) 



[-3± V9 + 6w], (78) 



in which we assume w ^ and also the condition (|65p . Putting eq. (|78p into 
eg .([75]). gives 

where 

.M w ) = -[3(w+ 1) T V / 3(2w + 3) ]. (80) 
w 

In this way the exact solution of the Brans-Dicke model is found as follows: 

H(z)=H(z )(l + z) 9 ^) . (81) 
The equation of state parameter u can be found as 

which is constant in time. So it is natural that it does not cross the phantom- 
divide-line, and lu(z) = —1 occurs when g^(w) = 0. <7-(w) = results in 
w = which is not acceptable, and <?+(w) = results in w = —4/3 which is in 
accordance with the perturbative result obtained in eq. (|H7| . 

It is interesting to note that our exact solution (f8TT) coincides with O'Hanlon 
and Tupper solution: [23] 



a(t) =ao ( — J , (83) 



4>(t) = <Ao ( ^ ) , (84) 



with 



q± 



w 



In t-space, (T55]) results in 



3(w+l) T V 3 ( 2w + 3 ) ' 
IT ^3(2w + 3) 
3w + 4 



(85) 
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which by noting that g T = l/q± (see eqs. (|8"0")) and {55}), the u> in becomes 
the same as one in eci.([52"j). 

The only remaining case is w = —3/2. To see the behavior of the Brans-Dicke 
model atw = —3/2, we consider the Friedmann equations ((72^) -(f7i |l . From the 
first equation (|72|) . one finds 



1 4> _ 3 ± V9 + 6w 
H d> w 



(87) 



from which 



4> = -2H<t> , 
4> = -2H<j) + AH 2 (t> . 

Putting these equations back into the second and third Friedmann equations 
((75)1 and (|74")l . which are not independent at w = —3/2, it can be easily seen 
that these equations trivially hold for arbitrary H(t) function. Therefore we 
exclude w = —3/2 case in which all the H(t) functions are solutions of the 
Friedmann equations. 



7 The quantum corrections 

As it was discussed earlier, the classical Lagrangian 

£ = - (v^) 2 ) (89) 

can induce interesting quantum effects. The renormalizability requirements en- 
forces £ in the above action to be £ = 1/6. In £ = 1/6, the model has several 
other interesting properties, such as the conformal invariance, holding the Ein- 
stein equivalence principle, etc. [17,22]. 

Calculating the effective action of this model at one-loop level, results in 
some extra terms in the trace of the energy-momentum tensor, which is trace- 
less classically. For this reason, this effect is called the trace/conformal anomaly. 
These extra terms are [19,22]. 

T = b{F + ^DR) + b'G + b"DR , (90) 

where T denotes the trace of the energy-momentum tensor, F is the square of 
the 4d Weyl tensor and G is the Gauss-Bonnet invariant 

6 ( 91 ) 

G = R 2 ~ 4i? M1/ i? M1/ + R^yapR^ 01 ^ ■ 

For N scalars, Ni/ 2 spinors, Ni vector fields, N 2 (=0 or 1) gravitons and Nhd 
higher derivative conformal scalars, 6, V and b" are given by 

N + 6N 1/2 + 127Vi + 611ZV 2 - 8N HD 

^ ' (92) 

N + llN 1/2 + 62jV! + 1411JV 2 - 28N HD K ' 

360(4^) 2 ' ~ 
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Using eq. (|9"0|) for the case of FRW metric, it can be shown that the contribution 
of conformal anomaly to energy density and pressure are [24] 

P A = —^{b'{6a 4 H 4 + \2a 2 H 2 ) 

+ (|& + b ")[ a 4 {~6HH - 1%H 2 H + 3H 2 ) + 6a 2 H 2 } ^ 
-26 + 66' -36"} , 

and 

PA = b'[ 6iJ 4 + 8H 2 H + \(AH 2 + 8H) } 

+ (tt& + b ")[ ~ 2 H - \2HH - \8H 2 H - 9H 2 + \{2H 2 + AH) ] (94) 
3 ar 

-26 + 66' - 36" 
3a 3 ' 

The subscript " A" stands for " anomaly" . To consider the quantum effects in 
gravitational phenomena, one can add the above pa and pa to the Friedmann 
equations. In this way the contribution of quantum effects oiioj = -1 crossing 
can be studied. 



7.1 The coupled-quintessence model 

We first study how much the adding of the quantum terms pa and pa can 
change the uu = — 1 crossing of the coupled quintessence model discussed in 
section 5. As this model is a theory with one scalar field, one has N = 1 and 
N 1/2 — Ni = N2 = N HD = 0. Eq.dlll) then results in 

b = 3 ^120^' < 95 > 

The Friedmann equations of coupled quintessence model are eas. (|36|) - (138)1 . 
in which we put pa and pa as the source terms. At zero-order, eqs. (|93j) and 
(fMf . using ipg]). are 

p A (0) = -b'(8h 4 + 24/io/ir + (96) 

and 

PA (0) = 6'(6^ + A8h h 1 - \ + 24^ 2 ), (97) 
respectively. Therefore, instead of eqs. ([3"9")) and (|4"U|) . now we have: 

4 + 12e/i o 0o0o - 6^(1 - i4>l) ~ 26'(8^ + 24/io/n + i|) = , (98) 

a o 

and 

16 

- 2£/io0o0o + 2C0o0o + (2$ - 1)$ = b'{-2h* + 2Ah h 1 - — + 2Ah 2 ) , (99) 

a 

respectively. Eq. (|4"Tj) is not changed. As it is clear from the above equations, 
the situation is very different from classical relations (f3"9"l) - (|4"Tj) . Here we do 
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not have the relations like eqs. (|4"2")l and (|4"3"| which describe 0o in terms of cf>o, 
and therefore the constraint like eq. (f4"T|) does not appear. In other words, the 
u> = — 1 crossing can happen for any £ values. In fact, eq. (|98[) determines h\ as 
follows: 

^ 486% ' (100) 

which is generally different from zero for any £ values, including £ = 1/6. In this 
way we obtain an interesting result: The coupled quintessence model classically 
crosses the uj = — 1 line for £ > 3/16, but because of the quantum effects, the 
£ = 1/6 < 3/16 is also allowed. It is interesting to note that the quantum 
correction terms (the last two terms in the numerator of eq. (|100p ) are much 
smaller than the classical terms, since 

hi ~ 4 « h 2 . (101) 



7.2 The Brans-Dicke model 

In this case, the quantum induced energy density and pressure and (j9"7|) 
must be inserted into eqs. ([58|) - ([62|) . in order to find the influence of quantum 
phenomena on uj = —1 transition of Brans-Dicke model. Eq. (|60[) then becomes 

'2 i o 

- w^ + 6/i o 0o + 6hU + 2b'(8ht + 2Ah Q h x + = 0. (102) 

In contrast to classical case in which eq. tjSU]) . with p = 0, specifies ((/)/ (/>)<> in 
terms of "w" in eq. ([M|) . from which w = —4/3 becomes necessary for achieving 
u = — 1, see eq. (joT|) , and finally results in h\ = in eq. ([TTj) . here the eq. (pTj|) . in 
its new form (1102p . results in a nonzero h%: 

w0§/^q - 6fco0 o - 6feg0 o - 86 y (2/ig + 3/aj) 
hl 48Vh ■ (1 ° 3) 

This shows that the quantum effects can produce the non-zero hi which indicates 
the existence of u = — 1 crossing. 

There is a last point which must be noted. Comparing the Brans-Dicke 
Lagrangian ([57]) with Lagrangian ([89]) . which pa and pa in eqs. (|96]) and ([97]) 
are derived from, may lead us to conclude that it is not reasonable to use pa and 
Pa as the quantum corrections of ([571) . This is because these two Lagrangians 
are different, so the corrections come from ((Si?)) can not be used for (|57p . The 
answer is that if we apply the following change of field variable to eq. (f57|) : 



1 ^2 



4w 

the Brans-Dicke action then becomes 

i r .„ , — ( i 



(104) 



5 = - / d 4 x^ I ^-^ 2 - (V0) 2 ) . (105) 



By choosing £ = — l/4w, eq. (|105[) becomes the same as the Lagrangian 
so its quantum computations can be used for Brans-Dicke action. Restricting 
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ourselves to £ = 1/6, results in w = —3/2, which is not in the region of our 
interest. So in the physical region w > —3/2, the quantum phenomena can not 
change the noncrossing behavior of equation of state parameter from oj = — 1 
line. 



7.3 Deceleration to acceleration transition 

To study the influences of quantum effects on a — transition, the quantum 
terms pa and pa must be added to Friedmann equations j9]), (fT0|) and (I14| . and 
solved for the coefficients of the expansion H{t) in eq. (|3"Tj) . Eq. (|3"3"|) does not 
change, but the coefficient H2 in eq. l|3~i|) . for b = —36' case, is replaced by 

gfl = 12g D (2y + 3F RRm t »» + 5^ + 72if °^ - 3H ° F ^ , , 

(106) 

--F + 126'(3i? 4 --)] t=0 . 
z a 

In above equation, ao is the scale factor a(t) at the transition time. It is seen 
that in b' — > limit, eq. (|106p reduces to classical relation |3^|) . as it is expected. 
For special cases coupled quintessence and Brans-Dicke models, Hi can be easily 
calculated from eq. (|106|) . 



8 Conclusion 

In this paper, we consider the generalized scalar-tensor models as a class of 
modified gravity theories to describe two important cosmological transitions, 
the u> = — 1 and a — transitions. For both cases, we solve the Friedmann 
equations of generalized ST models by perturbative expansion of the Hubble 
parameter around the transition points. The expansion parameters are found 
consistently, which proves the existence of these solutions. 

Two specific examples are studied in details. The first one is the coupled 
quintessence model, which the scalar field is coupled to gravity nonminimally. 
It is shown that for all £ > 3/16 cases, the model has a u) = — 1 transition, 
a fact which can be seen by numerical solving of the Friedmann equations for 
several ^-values. The £ = 1 results are reported. The second example is the 
Brans-Dicke model, which our perturbative method shows no transition, a fact 
that is verified by solving exactly the Friedmann equations of this model. It 
is seen that the equation of state parameter of this model is constant, which 
naturally never cross the ui = — 1 line. Our investigation shows that "w" must 
satisfy w > —3/2. 

Finally we consider the quantum terms coming from the one-loop calcula- 
tion of a scalar field, coupled non-minimally to gravity. It is shown that these 
quantum field theoretical terms can, in general, change the u> = — 1 and a = 
crossing behaviors of the models. For instance it is shown that the coupled 
quintessence model with £ = 1/6, which classically can not cross the u> = — 1 
line, crosses this line because of quantum effects. For Brans-Dicke model, this 
phenomenon has no effect on transition of models with w > —3/2. 
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